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Abstract 

We have numerically studied the behavior of one dimensional tunnel junction 

arrays when random background charges are included using the "orthodox" 

theory of single electron tunneling. Random background charge distributions 

are verified in both amplitude and density. The use of a uniform array as a 

transistor is discussed both with and without random background charges. 

An analytic expression for the gain near zero gate voltage in a uniform array 

with no background charges is derived. The gate modulation with background 

charges present is simulated. 
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I. INTRODUCTION 



Since the discovery of correlated single electron tunneling many electronic devices ex- 
ploiting this effect have been proposed [|T]||. A gate is commonly used to control the device: 
In single electron transistors to modulate the current, in traps to pull electrons on to an 
island or in pumps and turn-styles to move electrons in a specific manner. The response 
to the gate is periodic in the induced charge on the gate capacitance, with a period of one 
electron. These devices are therefore very sensitive to charges trapped in the vicinity of 
the device, where there is a significant capacitive coupling to the islands. Randomly dis- 
tributed "background charge" can change the threshold voltage for injecting charge in the 
array and the gate voltage where the coulomb blockade is lifted. The background charges 
may be static, causing a shift of the device operation point, or dynamic, causing noise in 
the device 0,|. 

The effect of background charges has been observed in experiments with small arrays. 
The influence of the background charge shows up as an offset in the current versus gate 
voltage curves ||. It has also been shown that the threshold voltage could be maximized 
by tuning gates to each island and thereby compensating for the background charge ||. In 
single electron traps it was shown that the discrepancy between experiments and theory 
was likely to be due to random background charges J7|. Theoretical studies of the effect 
of random background charges in one dimensional arrays have examined how the threshold 



voltage is affected [§|§ and noise properties when background charges are included flG 



Here we discuss simulations of the transport properties of one dimensional arrays of series 
coupled tunnel junctions. The arrays are characterized by nearest neighbor capacitance C 
and gate capacitance C g . The simulations include random background charges which are 
characterized by their density and amplitude. We examine the effect of static background 
charges on the transport properties, including the gate voltage dependence, for various ratios 
C/C 9 . 
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II. BASIC EQUATIONS 



The one dimensional array consists of iV — 1 small islands connected in series by small 
capacitance tunnel junctions. An equivalent circuit is shown in fig. [1]. A gate voltage, 
V 9t i, is coupled through the gate capacitance C g ^ to each island. A transport voltage is 
applied across the array by the external voltages Vj, and Vr. The tunnel junctions are 
characterized by the capacitance Cj and the resistance Ri. If the tunnel resistance satisfies 
Ri S> h/e 2 = 25.81281... kfl, the number of charges on each island becomes sharply defined 
and the quasi-classical "orthodox" theory of single electron tunneling applies ]TJ . The theory 
predicts that electrons tunnel one at a time with a tunneling rate given by 

1 AF 

i ~ e 2 Ri 1 - exp(AF/fc 6 T) ' ^ 

where AF is the difference in free energy of the initial and final state. We are interested 
in the change in free energy when one electron tunnels through one junction in the array, 
say from island k to island k ± 1. In this case the free energy difference is given by the 
electrostatic potential difference of those islands before and after the tunneling event 

AF — F' — F — e -[{ Vk±l + <p' k±1 ) - {fp h + (2) 

In the case considered here, all junctions and gate capacitances are equal, Ri = R, Ci = C 
and C 5i j = C g . Charge conservation gives a relation between the potential on the islands 
and the charges on the islands: 

(2 + cy cv -<p 2 = qi /c + v L 

-tpi-x + (2 + C g /C) Vl - cp i+1 = qi/C, 2 < i < N - 2 (3) 
-if N -2 + (2 + C g /C)ip N ^ = q N -i/C + V R . 

where = erii + go,j + en 9 is the charge on island i, consisting of the induced background 
charge on that island g 0i i an d the induced gate charge en g . The equations @ must be 
inverted to find the potential on the islands. This inversion can be done exactly [|1^]. The 
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result of this inversion can be expressed in matrix form: <p = Rg, where the components of 
the symmetric matrix R are 

cosh(A r — \i — j\)X — cosh(iV — i — j)\ 
ij = 2sinhAsinhiVA ^ ' 

and A = arch(l + C g /2C) is known as the inverse of the charge soliton length introduced in 
ref. ||11|| . The charge soliton length, A -1 , is the screening length giving the distance over which 



is 



the potential from one excess electron is screened by the free charge in the ground plane. 
With equations (HI)-©, all tunneling rates can be calculated. A Monte Carlo scheme [|ll| 
used for the simulation of the flow of charges. The random numbers used to generate the 
background charge configuration and for the monte Carlo calculation are chosen to make a 



good quasi random series [13] . All results are calculated for T = 0. 



III. EFFECT OF RANDOM BACKGROUND CHARGES 

Charges trapped in the substrate or in the tunnel barriers induce an excess charge on 
the islands, qi — > + q^. The background charge qoi, is between -e/2 and e/2 because any 
value outside this interval would be compensated for by one electron or hole tunneling on 
to the island. 

Starting from a random background charge configuration, under certain conditions tun- 
neling events will lower the free energy by adding or rearranging electrons or holes. Since 
the electrostatic energy is proportional to the gradient of the potential squared, lowering 
the energy means smoothing the potential. This possibility to relax the random potential 
variations depends on the length of the array, the soliton length, and on the actual configu- 
ration of the background charges. In fig. Q the potential before and after relaxation is shown 
for two different values of A. One can see by looking at fig. 0a that for short soliton length 
the potential is nearly unchanged by the relaxation process. However for long soliton length 
(fig. 0b) we can see how the relaxation reduces variation of the potential. The potentials 
plotted in figure 2 for different C/C g are normalized by the threshold voltage. The threshold 
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voltage (derived below) is the proper measure of the Coulomb blockade for transport in an 
array. 

The random background charge configuration can be described by two parameters: the 
number of islands which have an induced charge qoj 7^ 0, or the density of induced charges 
p E {0,N — 1}, and the amplitude A of the induced charges q 0<i £ {— Ae, +Ae} with 
< A < 0.5. In fig. |H the mean value of the "relaxed" potential variation along the array 
(std <p), averaged over 1000 configurations of random background charges, is plotted versus 
p and A. The four plots in fig. § correspond to four different values of the soliton length 
A -1 . Here we can see that the smoothing of the potential becomes more effective when the 
soliton length increases. Comparison of figures |3]c and |^d show that when A -1 > N the 
value of (std ip) is weakly dependent on A. Further we can see that when A -1 > N and 
p > N/2, A > 0.25, the potential variations of the array are weakly dependent on A and 
p. This weak dependence can also be seen in fig. |] where the threshold voltage is plotted 
as a function of the amplitude (right panels) and density (left panels) for various A. The 
threshold voltage changes almost linearly with the amplitude and density of the background 
charges when A = 1 (figs. |]a and |^c), but when A = 0.01 the dependence is very weak for 
A > 0.25 and/or p > N/2 (figs. [|b and |]d). Thus we see that as far as the electrostatics is 
concerned, disorder due to random background charges may be effectively reduced by single 
charge tunneling if the soliton length is larger than the array length and the density and 
amplitude of random background charges is large enough. 

IV. GATE EFFECTS IN ID ARRAYS 

When a gate voltage is capacitively coupled to the islands of the array the current 
voltage characteristics are strongly affected around the threshold. In the case when the gate 
is uniformly applied (V g ,i = V g and C g ^ = C g ) and there are no background charges, the 
threshold voltage changes with induced charge as shown in fig. (|J). For appropriate values of 
N and A an array with uniform coupled gate capacitance has ideal transistor characteristics. 
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Such an array would make a good switch because of the initial steep dependence of Vt on 
n g , and the wide region of n g with a weaker dependence. 

The sensitivity to gate voltage is greatest when iV is large but N < A -1 . In this case the 



threshold voltage is large at zero gate voltage ||12|| , and the gate drives the threshold almost 
to zero. The initial linear decrease of V t with n g can be examined analytically because during 
this decrease the gate voltage is not large enough to bind any charges inside the array. When 
no charges are localized in the array (electrons, holes or background charges) it was shown 



in ref. [y| that the threshold voltage is determined by the tunneling rates through the two 
end junctions. With eq. (§) we find the relevant potentials to calculate the tunneling rate 
through the end junction 



N-l 



R 1A (l + CV L /e)+n g J2^t,j 

3=1 



c 



N-l 

Ri,iW/e + n 9 ^R l5 
i=i 

^0 = ^0 = V L- 



(5) 

(6) 
(7) 



Here we have taken Vl > 0, Vr = and V g < so in this case a hole will tunnel in from the 
left. Making use of the relations 

sinh(iV- j)X 



R-i, 



N-l 
3=1 



sinh(A^A) ' 
sinh[|(^-l)A] 



(8) 



(9) 



2cosh(iVA/2)sinh(A/2)' 
We obtain the tunneling rate from eq. (0). Setting the tunneling rate to zero and solving 
for Vl we get the threshold voltage: 



V t 



c 



sinh[i(iV- 1)A] sinh ATA/2 

— 7 

sinh A^A - sinh(AT - 1)A sinh 2 A/2 cosh(A^ - |)A 



sinh(AT - 1)A 



(10) 



If n g is set to zero in this expression we recover the expression for the threshold voltage found 



in 



T2"f. This equation explains the very steep slope of the decrease in threshold voltage, for 



large A 1 and N. 



Eq. ( |TOD can be used to calculate the voltage gain of the uniform array. If the array is 
current biased very close to zero and it is assumed that the gate voltage drives the threshold 
voltage to zero, the voltage gain is given by 

_ dV out _ C g dV t _ smh[|C/V - 1)A] sinhiVA/2 
V dV g ~ e dn g cosh(iV - §)A 1 ' 

The gain is plotted in fig. || for different values of A. The gain approaches 2 when AiV > 3, 
which resembles the voltage gain in the uniform single electron transistor where one needs 
C g ^> C to get a maximum gain of 2. 

Unfortunately this gain behavior near zero gate voltage is destroyed by random back- 
ground charges. Melsen et. al. || have shown that the mean value (averaged over background 
charge configurations) of the threshold voltage was independent of the gate voltage. Our 
simulations show that it is still possible to modulate the threshold voltage with gate voltage 
when background charges are present, however the modulation changes. The large peak and 
high sensitivity around zero gate voltage disappears, and the largest modulation can occur 
near any value of n g . The maximum modulation (V^| max — V t \min) also becomes smaller 
than in the no-background-charge case. In fig. [F] the maximum modulation (averaged over 
1000 background charge configurations) is plotted as a function of number of junctions for 
different values of A. As in the no-background-charge case, the modulation becomes larger 
with increasing N and A -1 . Here again, we see weak dependence of the modulation on N 
when NX^ 1 > 1. The error bars plotted for several points give the standard deviation of 
the modulation, or the spread in modulation depending on the particular configuration of 
background charges. Thus, static random background charges cause random shift in the 
operation point in gate voltage, reduce the maximum modulation with gate voltage, and 
cause significant variation in the gate modulation of the array. 

We have simulated the effects of random background charges on the Coulomb blockade 
in one dimensional tunnel junction arrays. We find that when A" 1 > iV the static potential 
along the array is smooth and does not depend very strongly on the density or amplitude 
of background charges. We derive an analytical expression for the voltage gain of an array 

7 



transistor, valid when n g is small and no background charges present. However, by simulating 
the behavior of the threshold voltage as a function of a symmetrically coupled gate voltage 
with background charges in the array, we conclude that random background charges render 
arrays useless in transistor applications. 

We gratefully acknowledge financial support from the Swedish TFR/NFR and the Eu- 
ropean Union through MEL ARI project 22953 CHARGE. 
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FIGURES 



FIG. 1. The schematic of a ID array of small capacitance tunnel junctions. 

FIG. 2. The potential in a 50 junction array as a function of the island number for A = 1 and 
A = 0.01. The solid line is before the free energy has been minimized by single charge tunneling, 
and the dashed line is after minimization. 

FIG. 3. The mean value of the standard deviation of the potential in the array of 100 junctions 
for different values of A. 

FIG. 4. The threshold voltage in a 50 junction array as a function of the amplitude, A (left 
panels) and the density,/? (right panels) of the background charges. The different curves on the left 
correspond to p = 49, 30, 15 and 5. The different curves on the right correspond to A = 1, 0.6, 0.3 
and 0.1. 

FIG. 5. The effect of the gate on the threshold voltage for an array with no background charges 
present, N = 100. 

FIG. 6. The voltage gain near zero gate voltage of an array plotted veurses N for different 



values of A. The analytic expresion is given in equ. 11 valid for arbitrary N and A. 



FIG. 7. The averagy modulation of the threshold voltage plotted versus N for various A. The 
errorbars give the standard deviation calculated from 1000 random background charge configura- 
tions. 
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